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Abstract 

The nonlinear magneto-optical response of s- and d-w&ve superconductors 
is discussed. We carry out the symmetry analysis of the nonlinear magneto- 
optical susceptibility in the superconducting state. Due to the surface sensi- 
tivity of the nonlinear optical response for systems with bulk inversion sym- 
metry, we perform a group theoretical classification of the superconducting 
order parameter close to a surface. For the first time, the mixing of singlet and 
triplet pairing states induced by spin-orbit coupling is systematically taken 
into account. We show that the interference of singlet and triplet pairing 
states leads to an observable contribution of the nonlinear magneto-optical 
Kerr effect. This effect is not only sensitive to the anisotropy of the gap 
function but also to the symmetry itself. In view of the current discussion 
of the order parameter symmetry of High-T c superconductors, results for a 
tetragonal system with bulk singlet pairing for various pairing symmetries are 
discussed. 
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I. INTRODUCTION 



The investigation of the symmetry of the superconducting order parameters is currently 
one of the most exciting problems in the field of high-T c research. 1-9 Many experiments 
suggest that the order parameter is anisotropic. 1-5 Furthermore, the observation of it phase 
shifts in corner junctions 6 ' 8 ' 9 is consistent with a d x i_ y i symmetry of the gap function, al- 
though another experiment 7 favors a more conventional s-wave pairing. On the other hand, 
the occurrence of a finite tunneling current perpendicular to the planes 10 seems to be incom- 
patible with rf-wave symmetry in a simple picture. In view of this debate, it is important to 
develop alternative experimental techniques which are able to discriminate between the var- 
ious symmetries of the gap function. In particular, it is interesting to measure the symmetry 
of the gap function without the necessity of a tunneling contact, which always includes the 
problems of a residual magnetic field, trapped flux in the tunneling circuit, or singularities 
in the supercurrent flow at the corner 11 . 

On the other hand, optical second harmonic generation has widely been used as a probe 
for two dimensional optical, electronic and structural properties at various interfaces such 
as solid- vacuum, solid-gas, solid-solid, solid-liquid, or liquid-gas 12 ' 13 . Recently, the nonlinear 
magneto-optical Kerr-effect has become a new nonlinear optical method for the investiga- 
tion of low-dimensional magnetic structures 14-16 . This effect is interface sensitive within the 
electric dipole-approximation and directly probes, unlike methods based on Raman scatter- 
ing, low-energy magnetic excitations (~ meV) using optical interband transitions (~ eV). In 
contrast to linear optical techniques, the nonlinear magneto-optical Kerr effect yields an ex- 
cellent signal to noise ratio, since it is free of linear nonmagnetic background radiation. This 
advantage manifests itself both in the large magnetic intensity contrast upon magnetization 
inversion, 

I(M) - I(-M) 
I(M)+I(-M) 

of about 40 % for Fe surfaces 17,18 compared to a linear effect of typically 0.1 % and in a 

(2) 

giant nonlinear Kerr rotation <p K which, depending on the chosen Kerr-configuration, can 
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be tuned up to 90° nearly at will 19-21 . This corresponds to an enhancement by 2 - 3 orders 
of magnitude compared to the usual linear Kerr rotation, which is even reduced in thin 
films. Thus, it allows for the determination of the magnetic interface symmetry, including 
the magnetic "easy axis" 22 . 

It is well-known that the symmetry of the nonlinear optical susceptibility is strongly 
affected by a magnetization or an external applied magnetic field. 14 Thus it is of consider- 
able interest to extend the symmetry analysis of the nonlinear magneto-optical response to 
the superconducting state for different symmetries of the superconducting order parameter. 
Therefore, we propose in this paper a new theory for the nonlinear optical response of a su- 
perconductor in the presence of a magnetic field, which is able to distinguish between certain 
symmetries of the gap function and might stimulate corresponding experiments. In the fol- 
lowing, we present our theory and show that indeed a symmetry dependent contribution to 
the magneto-optical response without tunneling contact results in optical second harmonic 
generation. Although it is, due to the gauge invariance, impossible to measure the phase of 
the superconducting order parameter without tunneling contact, it is still possible to mea- 
sure its symmetry, not solely its magnitude. This results from the interference of different 
pairing amplitudes in the dipole matrix elements of the three transitions in nonlinear optics. 
Due to the surface sensitivity of nonlinear optics for systems with bulk inversion symmetry, 
one can take advantage of the broken inversion symmetry at the surface. This is of interest, 
because Cooper pairing, together with the always present spin-orbit coupling, is then no 
more purely singlet- or triplet-like. The interference of the singlet and triplet pairing states, 
which is linear in spin-orbit interaction, leads to the symmetry sensitive contribution of the 
nonlinear optical response for systems in an external magnetic field. Note that the possible 
importance of mixed singlet and triplet pairing states at the interface of a tunneling contact 
between a heavy fermion and a conventional superconductor has already been discussed by 
Fenton 23 . In order to obtain a more systematic insight into these phenomena, we perform 
the description of a superconductor at the surface of a bulk inversion symmetric system, 
including spin-orbit interaction and give a group-theoretical classification of the irreducible 
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representations of the gap function. This is of importance for a detailed calculation of the 
tensor elements of the nonlinear optical susceptibility which is performed in the second part 
of the theory. Finally, we present our results concerning the symmetry dependence of the 
corresponding experiment and discuss in detail the differences of the line shapes of the opti- 
cal spectra. We find that it is possible to discriminate between an isotropic s-wave, a s x 2 +y 2- 
or d x 2_ y 2- wave and a d xy - wave. However, no symmetry dependent differences between the 
s x 2 +y 2- and the d x 2_ y 2- waves occur, which are actually the most discussed symmetries of 
the high-T c systems. Nevertheless nonlinear magneto-optics is shown to be an alternative 
and complementary method to gain insight into the symmetry of the superconducting gap 
function. 

II. THEORY 

The strategy of this section is to calculate the effects of s- and d-wave superconductivity 
on the nonlinear optical susceptibility tensor x^ 2 \ which are due to (a) the modification 
of the bandstructure and thus the resonance denominators and (b) due to the symmetry 
of the superconducting order-parameter affecting the wave functions of the optical electric- 
dipole transition matrix elements. Thus, we proceed as follows: (i) First we set up the 
superconducting BCS-type Hamiltonian and perform its group theoretical analysis. In order 
to obtain the desired sensitivity of nonlinear optics to the symmetry of the gap function this 
requires, as it will turn out, the simultaneous absence of inversion symmetry, presence of an 
external magnetic field breaking time reversal, and spin-orbit interaction coupling singlet and 
triplet pairs, (ii) Making use of this symmetry classification and the mentioned constraints 
we then calculate the nonlinear magneto-optical response in s- and d- wave superconductors 
from the appropriate current-current-current correlation function and propose a suitable 
experimental geometry for the observation of this new nonlinear magneto-optical effect in 
superconductors. 

In our theory for the nonlinear magneto-optical response of unconventional superconduc- 
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tors, the superconducting state is described within the BCS theory 24 . The corresponding 

— * 

Hamiltonian with arbitrary pairing symmetry and in a magnetic field h is given by 

# = E*l^iA, (i) 

with the four-component Nambu spinor = (ck M |, Ck M j., cLwf c -kuj.)- Here, c kMcr is the 
creation operator of an electron with momentum k, band index /x and spin a. The (4 x 4) 
matrix 7i kM can be expressed in terms of (2 x 2) block matrices: 



(2) 



y -A*_ kM -e^a° + h ■ 3 j 

The block matrices are expanded in terms of the unit matrix a° and the vector of the 
Pauli matrices a. This notation is close to that of Sigrist and Ueda 25 . The symmetry of the 
superconducting order parameter is characterized by the gap function A CTO -/ ;kM = (c ko -^c_ kcr ' M ). 
We neglect any diamagnetic, i.e. Meissner effect of the magnetic field, but assume a large 
penetration depth at the surface and no influence of the vortex structure to the optical 
spectrum. This seems to be reasonable at least for the excitations in the interband regime 
discussed in this paper. A CT(T / ;kM is decomposed in the usual way in singlet states (A kM = 

— * —* 

A° k/ J and triplet states (<i kM = —d_ k ^): 

A kM = (A kM <r° + ■ 3) ia y . (3) 

— * 

The symmetry of A k ^ and with respect to the transition from k to — k is a direct 
consequence of the Pauli principle. Since we consider the states at a surface, k refers to the 
two dimensional in-plane momentum. 

Below the transition temperature T c , the symmetry of a system is reduced compared 
to the high temperature phase. The symmetry group G of the high temperature phase is 
determined by the symmetry operations which keep the Hamiltonian for A k = invariant. 

— * 

We consider a system which is for h = invariant with respect to the group 

G = g x K x U(l) , (4) 



where g, K and U(l) are the point group, time reversal operation and the gauge group 
of multiplication of electron creation operator by an arbitrary phase, respectively 25-30 . In 
the ordinary case the normal-state gauge symmetry is broken at the superconducting phase 
transition, i.e. the residual symmetry group is g x K. This is called a conventional super- 
conductor. In unconventional superconductors however, the symmetry is lower than g x K. 
At the transition temperature, the BCS - gap equation is an eigenvalue equation and con- 
sequently, an eigenvector belongs to one of the irreducible representations V of the 
group G. If V is the unit representation Ai (or Ai g for systems with inversion symmetry) 
conventional superconductivity occurs. In all other cases the superconductivity is uncon- 
ventional. In order to discuss the various symmetry states of the order parameter, one has 
to generate all irreducible representations of the gap function, where, due to spin-orbit cou- 
pling, the spin degrees of freedom cannot be transformed independently from the spatial 
(orbital) coordinates. For various point groups this symmetry classification has been per- 
formed 25-30 . In all these cases the inversion operation Cj is an element of the group G. Since 
we are interested in the investigation of superconducting properties with surface sensitive 
nonlinear optical experiments, we have to take the effect of broken inversion symmetry into 
account. In order to be specific, we consider the surface of a tetragonal system (bulk point 
group D 4 h) with residual point group C± v . This group has five irreducible representations: 
four of dimension one (Ai, A2, &i and £> 2 ) and one of dimension two {£). The isotropic s Q 
and and anisotropic s x 2 +y 2 s-waves transform as A\, the d x i_ y i -wave as B±, the d xy -wave as 
£> 2 , a d x 2_ y 2d xy -wave as A2 and the p x , p y waves as E. 

In order to classify the irreducible representations of the gap function, we have to analyze 
the transformation properties of A^^^. Applying an element R of the point group to A^, 
the following transformation of the singlet and triplet part results: 



Here, T> R is the representation of R G G which transforms the coordinates. If one considers 
the transformation of a Pauli spinor with respect to a combination R = R Q Ci of the inversion 




(5) 



6 



operation and a rotation R Q , only the rotational part has to be applied to the spinor, i.e. 

the representation of R in spin space is T>^ 2 \ Consequently, for the vector in spin space 

dk, the representation = T>^, where the inversion operation is replaced by the identical 

transformation, has to be applied. Therefore, one finds in the bulk of a system with inversion 

symmetry: C{ A kfl = — d kfl ■ tfj icjv, since the vector cT is not affected directly by the 

inversion operation. The minus sign results from the inversion of k to — k. Consequently, in 

the bulk , the singlet and triplet part belong to different irreducible representations and either 

singlet or triplet superconductivity occurs. In distinction to this, a coexistence of singlet 

and triplet pairing states is possible for systems without inversion symmetry, i.e. at the 

surface. Now, the in-plane inversion operation can be realized by a rotation (rotation by tc 

— * — * 

with z-axis as rotation axis). This rotation transforms the vector d to — d and the minus sign 
of the transformation k — > — k is canceled. Consequently, the irreducible representations of 
the gap function at the surface contain both, singlet and triplet parts. 

From these considerations one obtains the irreducible representations of the gap function 
from the simultaneous Clebsch-Gordon coupling of orbital and spin degrees of freedom. The 
results for the simultaneously occurring singlet and triplet part of the pairing amplitude 
are given in the table. Sigrist and Rice 30 calculated the irreducible representations of the 
tetragonal group D^. Since = C± v x Cj, it is straightforward to check the above 
results by reducing the subduced representations of D ih . One finds: A\ = A\ g © A 2u , 
•A2 = A 2g © Aiu, Bi = Big © B 2u , B 2 = B 2g © B iu and S = S g © S u which leads to the 
results of the table. Here, g and u refer to the irreducible representations of with even 
and odd parity. The irreducible representations of table I are the possible symmetry states 
of a superconductor on the surface of a tetragonal system and with spin orbit interaction. 
If one neglects the spin orbit coupling, the spin and orbital degrees of freedom transform 
separately and the singlet and triplet pairing states decouple again. Therefore, for a bulk 
singlet superconductor, the simultaneously occurring triplet part at the surface is of the 
order of the spin-orbit interaction, and vice versa. In heavy fermion systems the spin-orbit 
interaction is large, but even in transition metals one expects this quantity to be of the order 
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of 50 meV which, although being smaller, has nevertheless dramatic consequences such as 
the reorientation of the magnetic easy axis in thin ferromagnetic films upon the increase of 
the film thickness or the rise of the temperature. Finally, we expect the spin-orbit induced 
triplet part to be observable if one considers a surface sensitive experiment such as second 
harmonic generation. 

Based on these group theoretical classifications, we calculate now the nonlinear magneto- 
optical susceptibility tensor of a superconductor and focus on the interference of the simulta- 
neously occurring singlet and triplet part of the gap function. The optical response in second 
harmonic generation can be obtained from the nonlinear current-current-current correlation 
function: 

r°° de r°° de' de" , „ 

Xa/3 7 (q,^)=/ — / — / — ^a/3 7 (q,e,e ,e ) x 

J-oo 7T J-oo 7T J-oo 7T 

f{e")-f(e') _ /(e')-/(e) 
u>+i6—e" +e' u>+iS—e'+e 

2(u + i8) - e" + e 1 ' 

where /(e) is the Fermi function and the spectral function 7 a/ g 7 (q, e, e', e") is given by 

J a/ 3 7 (q, e, e', e") = Tr ( J^ a Q{e)J fi pQ{e')J tri Q{e")) • ( 7 ) 

J qa is the a-th component of the current operator 

J<i = X] ik^c ] j c+ a CTM c k _a (TI , , (8) 

and g(e) = — ^lm(e + i5 + H)^ 1 is the density of states matrix with Hamiltonian H of Eq. 2. 
The trace has to be performed with respect to all single particle states, i.e. the momentum 
(k), band (/x, is), spin (a) and Nambu degrees of freedom. In the following, we consider 
only interband transitions /i ^ and the limit of the dipole approximation q — > can be 
performed without special care for plasmonic excitations. 

The evaluation of the trace of Eq. 7 is straightforward with the knowledge of the unitary 
transformation U^, which diagonalizes H and g(e). In the following we discuss for a 
bulk singlet superconductor in the limit of weak spin-orbit interaction A s o . and weak external 
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field because the phase sensitive contributions of the optical response vanish if either A s . a 
or the magnetic field is zero. For weak external magnetic field, the eigenvalues are given by 



= ±h ± £, 



(9) 



where 



(10) 



h is the absolute mangnitude of h and tr denotes the trace in spin space. A kAt belongs to 
one of the irreducible representations of table I. Analogously, the unitary transformation is 
given by 



(11) 



where and are the transformations which diagonalize TL^ for h = and = 0, 
respectively. This is correct up to first order in X s o and h. The zero-field transformation 
can be expressed in terms of (2 x 2) matrices w kM and -0 kM 25 : 



( . . \ 



(12) 



where 



(13) 



and 



i>kn — — A kM / £ kM . 



(14) 



Similar to the gap function in Eq. 3, these (2 x 2) matrices are expanded in Pauli matrices, 
leading to a singlet part i> k and triplet part of v^. The transformation is determined 
by a rotation in spin space 



exp(-ia-a) 

exp (id ■ a* 



) 



(15) 
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with rotation axis e a = — (e 2 x eh)/\e z x e\\ and angle cos(2a) = e 2 • (a = ae a ). is 

— * 

the unit vector in the direction of h. Using this transformation, we can diagonalize the 
Hamiltonian leading to the quasiparticle spinor 



k/i — ^k^^k^ 



(16) 



Here the components of are the destruction operators of the eigenstates with eigenvalues 
given in Eq.( 9). 

In order to express the current operator in terms of the Nambu spinors, we have to 
consider the behavior of the matrix element jvp,v under the simultaneous transformation 
k — > — k and (/i, v) — > (z/, p): 



jiv r 



(17) 



Although the phase factors depend on the choice of the phase of the Wannier functions, 
all observable quantities like Xa/3-y(<l, w) are independent of this choice. Now, the current 
operator can be expressed in terms of the Nambu spinors: 





Hi/ ±.o 



(18) 



V -j£d° I 



and we find for the nonlinear spectral function of Eq. ( 7) : 



1 



Iaf3y(t,£ , C ) — „ X/ ^k^i/iki//tik/t/i 

\ifJLVK, 

xTr' {AC ei{e) Ml K eS(0 AC e£(0> • 



(19) 



Here, AC and f?k(e) are (4 x 4) matrices in Nambu and spin space, whereby 

Uvu (20) 







V 



-p^a° 

results from the transformation of Eq.( 18) into the quasiparticle representation, where f?k(e) 
is diagonal. Consequently, the trace Tr' has to be performed with respect to the Nambu and 
spin degrees of freedom. 
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Nonlinear optics and in particular nonlinear magneto-optics offers a unique method to 
probe low-lying excitations close to the Fermi-level with optical photons, since SHG, which 
involves three photons, takes advantage of an additional degree of freedom that is absent 
in linear optics. Thus it allows to use conventional monochromatic, intense, and tunable 
pulse laser sources in the ps to fs regime such as mode-locked Ti-sapphire lasers. Besides, 
in nonlinear interband optics, there is no collective plasmon background which is material 
insensitive and may severely hamper the interpretation of linear optical experiments. 

Therefore, we restrict ourselves to a special interband excitation process. We consider 
the transition from the initial state % with energy Ei ps — 3eV, below the Fermi energy to the 
intermediate state s at the Fermi level (which is the only superconducting state) and to the 
final state / with energy Ef ps 3eV above E F , i.e. we consider /i — f, k — s, and v — %. A 
possible, but not necessary, origin of the states % and / might be due to the Mott-Hubbard 
splitting of the hybridized Cu 3d x 2_ ?/ 2 and O 2p x ( J/ ) orbitals. Since the intermediate state is 
the only state with superconducting coherence, we skip the band index of the matrices 
and -Ok- Performing finally the traces in Eq.( 19), we obtain: 

= + X$» + 0(h 2 ) , (21) 

where 

X$» = E J^jLjlf (K| 2 Gi(k,c) + K| 2 G 2 (k,c)) (22) 

k 

is the zero field susceptibility tensor in second harmonic generation within the super- 
conducting state which gives already a contribution without spin-orbit interaction. The 
nonvanishing tensor elements for Xa^i^) are the same as m the normal state G 
{zzz, zxx, zyy, xzx, xxz, yzy, yyz}). This results here from the transformational properties 
of the three matrix elements jkfiJkisJksfi which transform as the corresponding combination 
of the coordinates x Q xpx 1 , even if a single matrix element e.g. does not transform like 
x a . The functions Gi(2)(k, u) result from the numerous combinations of Fermi functions 
and energy denominators which occur by performing the traces in spin and Nambu space. 
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Xa^y( u ) * s the new contribution of the magneto-optical Kerr effect in the superconducting 
state. This can be seen from the symmetry relations of the magneto-optical susceptibility 

xSrH = E Jifi JL jLf KY v k ■ 4 F(k, u) , (23) 

k 

which depends on the superconducting gap function not only through its magnitude but 
also through A£ itselfe. Due to the additional triplet part however, the result is still gauge 
invariant. The function F(k, u) correspond to the Gi( 2 )(k, u) for X*£py( u )- Considering 
a magnetic field parallel to the x-axis (in plane), the nonvanishing elements of Xap-y( u ) 
are: a(5^ G {yyy,xxy,xyx,yxx, zzy, zyz,yzz}. This results from the combination of the 
transformation properties of the normal state matrix elements and of the symmetry sensitive 
term (v£)*Vk • 4- Using Eq.( 14), the k-dependence of the latter results from the irreducible 
representations given in table I. 

The above matrix elements lead to a rotation of the polarization of the incident light due 
to the interference of the singlet and triplet states at the surface of a superconductor. Thus 
nonlinear magneto-optics, unlike linear optical probes, provides indeed an optical method 
to discriminate different superconducting pairing symmetries by exclusively employing the 
effect of optical photons to low energy excitations. In the next paragraph, we discuss the 
numerical results of the above model bandstructure for the most realistic experimental setups 
and show how experiments can distinguish between certain symmetries of the gap function 
using nonlinear magneto-optics. 

III. RESULTS 

In this section we discuss the numerical results obtained for the tensor elements xf) z °f 
second harmonic generation without magnetic field and Xy h J z which gives rise to the rotation 
of the polarization plane for an applied magnetic field parallel to the x-axis. For simplicity we 
neglect the dispersion of the initial and the final state and consider solely the k-dependence 
which results from the superconducting gap function. The momentum summations are 
performed within the two dimensional Brillouin zone using 81 x 81 k-points. 
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The calculations are performed for a magnetic field of 9T (corresponding to a field 
induced band splitting of 0.5 meV) and a temperature of 1.5 K. The magnitude of the 
singlet part of the gap function is assumed to be 5meV. Furthermore, the magnitude of 
the dipole matrix elements is estimated to be 10 -11 m. All results presented here are not 
sensitive to the specific set of parameters chosen, but are typical for reasonable values of the 
corresponding energy scales. This was checked by systematically varying the dependence 
of the nonlinear susceptibility on the magnitude of the gap function, the magnetic field, 
the position of the initial and final states Ei and Ef, the temperature and the linewidth 
broadening S. 

In Figs. 1(a) and (b), we show w 2 Imx^(w) and oj 2 lmx^ z {oj) for an isotropic s-wave. 
For the conventional SHG, we find a line shape similar to the real part of a Lorentzian, 
which is typical for a three level system discussed in this paper. More interestingly, the 
fine structure of the peak, shown in the inset of Fig. 1(a) clearly shows the energy scale 
of the superconducting gap. Comparing this behavior with the Kerr signal w 2 Imx^ 2 (w) of 
Fig. 1(b), one finds that the interference of the singlet and triplet pairing states leads to a 
line shape with several pronounced zeros and with a fine structure that yields, besides the 
energy scale of the superconducting gap, also excitations which result from the magnetic field 
splitting. In all our calculations, this line shape was exclusively observed for an isotropic 
s-wave and can be considered as a fingerprint of this symmetry. 

In Figs. 2 and 3, the corresponding results for the anisotropic s-wave and the d x 2_ y 2-wave 
are shown. Although the result for the conventional SHG is similar to that of the isotropic 
s-wave, a totally different line shape of the Kerr signal results. This is due to the symmetry 
dependent prefactor in Xyzz( u ) an< ^ can be used to discriminate these two symmetries from 
the isotropic s-wave. Furthermore the fine-structure of these two symmetries is very different 
from the isotropic s-wave. Due to the occurrence of nodes in the gap, not only a peak but a 
whole broad band between 3eV and 3.01 eV is observable. This range is surprisingly given 
by twice the superconducting gap magnitude. Unfortunately, there are only slight differences 
between the two symmetries shown in Fig. 2 and 3. This is due to the similar k-dependence 

13 



of the triplet part given in table I. Only the fine structure of the peaks displayed in the 
insets of Fig. 2(b) and 3(b) exhibits a clear difference, where the anisotropic s-wave has a 
clear zero at 3.01 eV which is more or less smeared out for the d x 2_ y 2 wave. 

In this context it is of importance to compare these results to the usual nonlinear Kerr 
effect in the normal state in an external applied field which can be easily estimated using 
previous results by Pustogowa et al. 31 as 

~)Cyzz, normal ^ ^s.o 

Here, the incident photon energy Tiuj at resonance is 3 eV, while the energy splitting J{h) 
caused by the external magnetic field h is 0.5 meV (see above). For this small value of 
J{h), the formula given in appendix C of the paper by Pustogowa et al. yields a linear 
dependence. Thus, we find yzz, " ° rmai ~ 2.7 • 10 6 due to the absence of a spontaneous 

\yzz 

magnetization. From this estimate, it follows that the observability of the new contribution 
to the nonlinear Kerr effect is clearly guaranteed for the anisotropic s-wave. Although, 
the intensities of the other symmetries are smaller than the estimated value of the usual 
nonlinear Kerr effect in both the normal and superconducting state, we believe that this 
effect is still observable for the following reason: Due to the neglect of the dispersion of the 
states in our model bandstructure, the disapperance of the signal results from cancellations 
of contributions of the order of magnitude of the anisotropic s-wave in the k-summation. 
This is due to the artificially high symmetry of this bandstructure. The simplified description 
of the high-T c materials is used because the aim of this paper is to demonstrate the strong 
interdependence of the lineshape of the nonlinear magneto-optical susceptibility and the 
symmetry of the superconducting state. A more realistic bandstructure immediately leads 
to larger intensities of x y h ) z i while keeping the characteristics of the line shapes of the spectra. 

In Fig. 4 we finally show our results for the d xy symmetry. Although, this symmetry 
does not seem to be the most probable candidate for the high-T c materials, it shows most 
clearly the symmetry dependence of the nonlinear magneto-optical Kerr effect. In contrast 
to the anisotropic s-wave and the d x 2_ y 2 wave, there occurs a sign change between the 
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magnetic and nonmagnetic optical spectrum. This is observable since the sign of the Kerr 
spectrum determines the direction of the rotation of the polarization axis, i.e. the Kerr 
angle. Furthermore, the satellites shown in the inset of Fig. 4(a) and (b) cover only the 
range from 3eV to 3.005 eV, i.e. only one times the gap magnitude. 

For the existence of a finite Kerr signal, it is necessary to break time reversal symmetry 
and to apply an external magnetic field. This enables one to keep any direction of the field 
fixed and to study the anisotropy of the effects discussed in this paper. However, due to the 
strong but short-ranged antiferromagnetic correlations it might also be possible to take ad- 
vantage of the locally broken time-reversal symmetry of the high-T c materials. Since a finite 
Kerr signal is expected for certain long-range ordered antiferromagnets 32 , a pump-and-probe 
experiment (on a time scale faster than the average lifetime of the local spin configurations 
T S pi n ~ 10 -12 — 10~ 13 s) could be able to resolve the influence of the neighboring spins on 
the site which is excited by the optical excitation. Furthermore, for a practical realization 
of the experiment, one has to take into account the possible heating effects of the incoming 
light on the sample, which might lead to a local disappearance of the superconducting state. 
In view of the comparable excitation energies of the nonlinear magneto-optical Kerr effect 
in ferromagnets, we believe that this effect is of minor importance, since for fs laser pulses 
sample heating is of the order of 10 K and thus negligible. This simple estimate is readily 
obtained from the comparison of laser heating with ns pulses 33 yielding intensities of 100 
MWm~ 2 and temperature rises of the order of 100 K and typical fs measurements of the 
nonlinear Kerr-effect 21 which operate at laser powers as low as 100 mW focused on spot 
diameters of 100 /xm. 

In conclusion, we presented a theory for the nonlinear magneto-optical response of su- 
perconductors. The surface sensitivity of this experiment is of particular interest for the 
simultaneous occurrence of singlet and triplet pairing amplitudes. Therefore, a suitable ma- 
terial for this experiment is the Bi 2 Sr 2 CaCu208 compound, where almost no reconstruction 
of the cleaved surface (Bi-0 layer) occurs, and where the existence of the superconducting 
state in the upper CuC>2 layer was clearly shown in photoemission experiments. 2 . Since, due 
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to difficulties in the manufacture of the tunneling contacts, all corner junction experiments 
where so far performed with YBa 2 Cu 3 7 _5 samples, this gives also insight into the symmetry 
of another class of cuprate compounds. Note that furthermore the surface sensitivity does 
not depend on the actual depth beneath the surface where superconductivity starts, since 
the electronic symmetry of the superconducting state surface is of relevance in this context 
which may or may not be perfectly identical to the physical surface. On the other hand the 
new nonlinear magneto-optical effect proposed in this paper may also be of considerable im- 
portance for interfaces between s- and rf-wave superconductors. Furthermore, we performed 
a group-theoretical classification of the irreducible representations of the gap function at the 
surface of a tetragonal system. Based on this theory, we showed that the interference of sin- 
glet and triplet pairing states in a magnetic field leads to a new contribution to the nonlinear 
magneto-optical Kerr signal which is sensitive to certain symmetries of the superconducting 
order parameter rather than only to its magnitude. This enables us to give the basic line 
shapes of the corresponding optical spectra and to show that it is possible to discriminate 
an isotropic s-wave and a d xy wave from the anisotropic s x 2 +y 2 and d x i_ y i waves. Unfor- 
tunately, it is not possible to discriminate between the two latter symmetries which seem 
to be the most probable symmetries of the high-T c materials. Nevertheless, we believe that 
the nonlinear optic can yield information complementary to the tunneling experiments and 
might be also of importance in view of the application on heavy fermion superconductors, 
where it was also manifested that the superconducting state is anomalous 34-36 . 
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TABLES 



irred. representation 


singlet part 


triplet part 


Ai 


const. , cos(k x ) + cos(ky) 


e x sm(k y ) - e y sm(k x ) 


A 2 


(cos(k x ) — cos(ky)) sin(k x ) sin(k y ) 


e x s\n(k x ) + e J/ sin(fcj / ) 


Bi 


cos(k x ) — cos(k y ) 


e x sm(k y ) + e y sin(k x ) 


B 2 


s'm(k x ) sin(k y ) 


e x s\n(k x ) - e y sm(k y ) 


£ 


sm(k x ) sin(fc 2 ), sin(fe^) sin(A; 2 ) 


e z sm(k x ), e z sin(k y ) 



TABLE I. Singlet and triplet parts of the gap function for the irreducible representations of 
the point group C4. k x and k y are the components of the in plane momentum. k z represents an 
additional quantum number, which changes sign if one interchanges the two paired electrons and 
which is related to the layer index at the surface. 
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FIGURES 

FIG. 1. w 2 Imxi°z(w) (a) and uj 2 lm.x y h )z{u) (b) for an isotropic s-wave symmetry of the gap 
function. Note the very different line shapes of the two experiments, typical for the isotropic 
s-wave. The insets show the fine structure of the results near the excitation energy of 3 eV. For 
better visibility the spectra are artificially broadened by a Lorentzian width 5 of 50 meV in the 
main figures and of 0.5 meV in the insets throughout. Consequently, the inset peak heights differ 
from those of the main figures. 

FIG. 2. As in Fig. 1, but for an anisotropic s x 2 +y 2 wave. Note the difference of the magnetic 
line shape compared to the isotropic s-wave and the broad fine structure up to two times the 
magnitude of the gap, shown in the inset. 

FIG. 3. As in Fig. 1, but for a d x 2_ y 2 wave. The spectra are very similar to the s x 2 +y 2 wave, 
which results from the similar k-dependence of the triplet pairing amplitude, given in the table. 
Slight differences occur for the fine structure near 3eV. 

FIG. 4. As in Fig. 1, but for a d xy wave. In distinction to the s x 2 +y 2 and d x i_ y 2 wave, the 
signs of the nonmagnetic and magnetic spectra are different. Furthermore the satellites shown in 
the inset occur only up to the magnitude of the gap function. 
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